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Abstract. We review the main theoretical aspects of the structure formation 
paradigm which impinge upon wide angle surveys: the early universe generation of 
gravitational metric fluctuations from quantum noise in scalar inflaton fields; the 
well understood and computed linear regime of CMB anisotropy and large scale 
structure (LSS) generation; the weakly nonlinear regime, where higher order per- 
turbation theory works well, and where the cosmic web picture operates, describing 
an interconnected LSS of clusters bridged by filaments, with membranes as the in- 
trafilament webbing. Current CMB+LSS data favour the simplest inflation-based 
ACDM models, with a primordial spectral index within about 5% of scale invariant 
and Q.a ~ 2/3, similar to that inferred from SNIa observations, and with open CDM 
models strongly disfavoured. The attack on the nonlinear regime with a variety of 
N-body and gas codes is described, as are the excursion set and peak-patch semi- 
analytic approaches to object collapse. The ingredients are mixed together in an 
illustrative gasdynamical simulation of dense supercluster formation. 

1 Introduction 

By their very nature, wide angle surveys, whether of CMB anisotropics or 
galaxy redshifts, are designed to probe long wavelengths in the universe, i.e., 
low comoving Fourier wavenumbers k (Fig. |l]). We are hoping for simplicity 
of interpretation, best if the observables are probes of linear physics, harder 
to disentangle if the probe is of the dissipative chaotic physics characterizing 
the highly nonlinear regime, and hopefully tractable in the mildly nonlinear 
regime. The goal is to extract the underlying linear gravitational metric fluc- 
tuation field $ which characterizes the initial conditions and the cosmic pa- 
rameters that define its evolution, whatever the specific CMB or density-based 
observable. 

A theorist would argue that it is best to go for CMB probes, using spectac- 
ular forecasts J|, H) of the percent level accuracy that cosmic parameters can be 
determined to because of the sensitive dependence of the anisotropies on the 
linear gas physics at photon decoupling at z ~ 1000. However, there are sig- 
nificant near-degeneracies among cosmological parameters such as A and the 
curvature, which require other experiments, often wide-angle, to break. Fur- 
ther, to realize the promise will depend upon how well foreground signals from 
the Milky Way and secondary anisotropy signals from nonlinear phenomena 
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Figure 1: This shows the bands in comoving wavenumber k probed by CMB 
primary and secondary anisotropy experiments and LSS observations. The 
nonlinear wavenumber fcjvi(t), where the mass density power crosses unity, 
propagates from high to low k with redshift leaving in its wake collapsed halos 
and at its bow the interconnected cosmic web pattern replayed with time on 
progressively increasing scales. The objects form from waves concentrated in 
the fc-space bands that their labels cover: the first star forming tiny dwarf 
galaxies appear ("1st *"), typically at redshift about 20; dwarf galaxies (dG); 
normal galaxies (gal); groups (gps); and clusters (els), where we are now. 
Equivalent mass scales are given above them. Secondary anisotropies arise 
only once matter has gone nonlinear. The range extends to the LSS regime 
because of clustering of the nonlinear objects, e.g., of clusters for the SZ effect. 

such as the Sunyaev-Zeldovich effect from high pressure gas and emission from 
dusty high redshift galaxies can be disentangled from our goal, the primary 
anisotropies. CMB surveys can deliver useful information all the way from the 
ultralarge scale structure (ULSS) realm "beyond" (our Hubble patch) through 
the highly linear VLSS regime well-probed by COBE, into the LSS regime. We 
are now getting increasingly good LSS information from current small angle 
CMB anisotropy experiments (bandpower estimates shown in Fig. [?]), we are 
about to take a significant step forward with balloon-borne and interferometer 
experiments, and a leap forward with the satellites MAP and Planck (lower 
right panel of Fig. ||) . 

Galaxy redshift surveys probe number density fields, so they are largely 
limited to probing the LSS band since the level of mass density fluctuations is 
apparently rather tiny in the VLSS band. The further complication inherent 
in biasing relative to the mass density field may still obscure the derivation 
of the gravitational potential field. In Fig. ||, this unknown is absorbed in a 
uniform-amplifier bias factor b g which the simplest biasing theories predict in 
the linear regime, but which we know could be quite environment and scale 



dependent in a way which may make reliable calculations difficult. 

We do think that for clusters the biasing has a simple physical explanation 
involving collapsed halos that allows for reliable calculations and extraction 
of the gravitational potential field from cluster surveys constructed using X- 
ray, optical or SZ data. It is this same simplicity of interpretation which has 
allowed the estimates of the mass power spectrum from cluster abundances 
(lower panels of Fig. ||). 

It seems that practical limitations will limit weak gravitational lensing and 
large scale streaming velocities to the LSS regime, even though they are direct 
probes of $ gradients. 

The length scales below LSS, labelled "gastrophysics" , are ones over which 
energy injection and propagation, through ionization and shock fronts from 
galaxies, quasars, etc. are expected to have had important effects. That 
is, even if the dark matter density power in the fc-bands in question (e.g., 
the galaxy band of Fig. [I] at z ~ 4) is linear, the gas effects can lead to 
severe and complex biasing or antibiasing over such scales. Thus, although 
probes based on wide field catalogues of high redshift galaxies and quasars, 
and on quasar absorption lines from the intergalactic medium, represent a very 
exciting observational frontier, it will be difficult for theoretical conclusions 
about the early universe and the underlying fluctuations to be divorced from 
these gastrophysical complications. 

1.1 Cosmic Structure and the Nonlinear Wavenumber 

Various (linear) density power spectra, V p (k) cx k i V^(k), are shown in Fig. |^. 
Many people plot V p /k 3 ~ kV<& instead. Here V$(k) = da%/dhxk is the 
rms power in each d In k band. In hierarchical structure formation models 
such as those considered here, the nonlinear wavenumber ki\rL(t), defined by 
Pp(kNL) = 1, where cr 2 p (k) = J V p (k)d\n.k, grows as the universe expands. 
kNh(t) was in the galaxy band at redshift 3 and is currently in the cluster 
band. 

At a given time t, there is a band in k extending just above fcjvx(i) which 
is primarily responsible for the nonlinear collapsed dark matter halos in the 
medium in hierarchical theories. Smaller halos that formed earlier from much 
higher k bands would have largely merged into the halos of relevance at epoch 
t, in a sequence of characteristic objects shown in Fig. |l|. If we denote the 
rms linear density fluctuation level for waves longer than k by a p (k), then the 
characteristic scale fc* which gives the peak of the Press-Schechter distribu- 
tion of collapsed objects, dH, co u(M) / d\n M , occurs at cr p (fc*) = f c , where f c 
is a suitable collapse threshold, the famous 1.686 for Q = 1. Associated with 
fc* is a characteristic mass^] M*. On the other hand, for mass scales below 
°> = /c/3 ~ 0.6 not even 3cr peaks will have collapsed at this resolution. As 

1 The relation between k and M = (4tt / '3) po I&, H is that the rms a p is the same for 
"top-hat filtering" on scale Rth and "sharp k-spacc filtering" on scale k~ . 
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Figure 2: These panels show how the shape of the linear density power spec- 
trum in the LSS regime changes as the cosmological model is varied. These 
are compared with the reconstructed linear power spectrum estimated by 
Peacockjl). It is unclear how seriously the discrepencies above 0.1( h _1 Mpc) _1 
should be taken. The two lower panels show how strong the overlap will be 
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between CMB and redshift surveys. Right focuses on V p , left on . A 
forecast of error bars for the number density power spectrum derived from 
the SDSS by Vogeley and Szalay, and current mass density power spectrum 
estimates from cluster abundances (CL) and large scale streaming veloci- 
ties (V), are contrasted with the dotted fc-space filters for MAP and Planck. 
Planck in particular encompasses well the two important scales which define 
the effective acoustic peak range for primary anistropies (those involving lin- 
ear fluctuations): the sound crossing distance c s r 7( ; ec at photon decoupling 
around redshift 1000; the width Ar 7dec of the region over which this decou- 
pling occurs, which is about a factor of 10 smaller, and below which primary 
CMB anisotropics are damped. Estimates of the linear $ power from cur- 
rent galaxy clustering data by Q] are compared with some sample (linear) 
COBE-normalized gravitational potential power spectra on the right. A bias 
is "allowed" to (uniformly) raise the shapes to match the observations. 



we describe later, this is the weakly nonlinear domain of the cosmic web of 
interconnected structure: the view of the structure smoothed on these scales 
reveals the characteristic patterns of filaments connecting clusters, and mem- 
branes connecting filaments. Voids are rare density minima which have opened 
up by gravitational dynamics and merged, opposite to the equally abundant 
rare density maxima, the clusters, in which the space collapses by factors of 
5-10 and more. 

At k > kpfi,(t), nonlinearities and complications associated with dissipative 
gas processes can obscure the direct connection to the early universe physics. 
Indeed, the reason galaxies are still around today when Unl is cluster scale 
is that, although the outer halos of the galaxies will have largely merged, gas 
cooling allows the baryons and some of the dark matter to concentrate and 
survive as independent beings in groups and clusters, breaking the hierarchy. 

Most easily interpretable are observables probing the linear regime now, 
k < fcjvx(to)- CMB anisotropies arising from the linear regime are termed 
primary; as Fig. [j] shows, these probe 3 decades in wavenumber, with the 
range defined by CMB damping rather than by any nonlinear effects at z ~ 
1000. LSS observations at low redshift probe a smaller, but overlapping, range. 
We have hope that z ~ 3 LSS observations, when fcjvx(i) was larger, can 
extend the range, but gasdynamics can modify the relation between observable 
and power spectrum in complex ways. Secondary anisotropies of the CMB, 
those associated with nonlinear phenomena, also probe smaller scales and the 
"gastrophysical" realm. 

1.2 Early Universe Connection and the Freedom in Inflation 

The source of fluctuations to input into the cosmic structure formation problem 
is likely to be found in early universe physics. The contenders for generation 
mechanism are (1) "zero point" quantum noise in scalar and tensor fields that 
must be there in the early universe if quantum mechanics is applicable and (2) 
topological defects which may arise in the inevitable phase transitions expected 
in the early universe. 

From CMB and LSS observations we hope to learn: the statistics of the 
fluctuations, whether Gaussian or non-Gaussian; the mode, whether adiabatic 
or isocurvature scalar perturbations, and whether there is a significant com- 
ponent in gravitational wave tensor perturbations; the power spectra for these 
modes, V$(k) ,Vi S {k) ,Vgw (k) as a function of comoving wavenumber k. As 
the Universe evolves the initial shape of V$ (nearly flat or scale invariant) 
is modified by characteristic scales imprinted on it that reflect the values of 
cosmological parameters such as the energy densities of baryons, cold and hot 
dark matter, in the vacuum (cosmological constant), and in curvature. Often 
observables can be expressed as weighted integrals over k of the final state 
power specta and thus can probe both Q's and initial power amplitudes and 
tilts. 

Many variants of the basic inflation theme have been proposed, sometimes 



with radically different consequences for V<s>{k) ~ fc 1- " 3 ^, and thus for LSS, 
which is used in fact to highly constrain the more baroque models. The most 
likely inflation possibilities are the simplest: (1) adiabatic curvature fluctua- 
tions with nearly uniform scalar tilt over the observable range, slightly more 
power to large scales (0.8 < n s < 1) than "scale invariance" [n s = 1) gives, 
a predictable nonzero gravity wave contribution with tilt similar to the scalar 
one, and tiny mean curvature (Qtot ~ 1); (2) same as (1), but with a tiny 
gravity wave contribution. 

An arguable rank-ordering of the more baroque inflation possibilities (see 
e.g., jj], [| for references) is: (3) same as (1) but with a subdominant isocur- 
vature component of nearly scale invariant tilt (the case in which isocurvature 
dominates is ruled out, but see Figs. ||,||): (4) radically broken scale invariance 
with weak to moderate features (ramps, mountains, valleys) in the fluctuation 
spectrum (strong ones are largely ruled out); (5) radical breaking with non- 
Gaussian features as well; (6) "open" inflation, with quantum tunneling or 
possibly a Hawking- Turok instanton resulting in a negatively-curved (hyper- 
bolic) space which inflates, but not so much as to flatten the mean curvature; 
(7) quantum creation of compact hyperbolic space from "nothing" . ft is quite 
debatable which of the cases beyond (2) are more or less plausible, with some 
claims that (4) is supersymmetry-inspired, others that (6) is not as improbable 
as it sounds. 

1.3 Cosmic Parameters and the CMB 

Even simple Gaussian inflation-generated fluctuations for structure formation 
have a large number of early universe parameters we would wish to determine: 
power spectrum amplitudes at some normalization wavenumber k n for the 
modes present, {V<p(k n ), Vi S (k n ), VGw(k n )}; shape functions for the "tilts" 
•{V s (fc), Vi s (k), v t (k)}, usually chosen to be constant or with a logarithmic cor- 
rection, e.g., v s {k n ),dv s {k n )/d\nk. (The scalar tilt for adiabatic fluctuations, 
v s (k) = dlnV^/dhik, is related to the usual index, n s , by v s = n s — 1.) The 
transport problem is dependent upon physical processes, and hence on physi- 
cal parameters. A partial list includes the Hubble parameter h, various mean 
energy densities {O, to t, ^b, ^a, ^cdm, ^/idm}h 2 , and parameters characterizing 
the ionization history of the Universe, e.g., the Compton optical depth tc 
from a reheating redshift z re h to the present. Instead of Sltot, it is becoming 
conventional to use the curvature energy parameter, = 1 — Qtot, thus zero 
for the flat case. In this space, the Hubble parameter, h = (X^(^j'-k 2 )) i 
and the age of the Universe, to, are functions of the fl/h 2 . The density in 
nonrelativistic (clustering) particles is Q n r = &b + Qcdm + Qhdm, denoted as 
well by Q m . The density in relativistic particles, J7 er , includes photons, rela- 
tivistic neutrinos and decaying particle products, if any. Q er , the abundance 
of primordial helium, etc. should also be considered as parameters to be de- 
termined. The count is thus at least 17, and many more if we do not restrict 
the shape of V$(k) through theoretical considerations of what is "likely" in 



inflation models. For example, the ratio of gravitational wave power to scalar 
adiabatic power is Vgw/V® « —(100/9)^/(1 — ft/2), with small corrections 
depending upon v s — Vt PJ. If such a relationship is assumed, the parameter 
count is lowered by one, and other restrictions on most likely behaviour or 
prior probabilities reflecting other observations or prejudices are used to re- 
duce the set further; e.g., the forecasts for the MAP and Planck satellites in 
[| used a 9 parameter set. 

For a given model, the early universe V§, is uniquely related to late-time 
power spectrum measures of relevance for the CMB, such as the quadrupolc 
or averages over ^-bands, and to LSS measures, such as the rms density fluc- 
tuation level on the 8h _1 Mpc (cluster) scale, ag, so any of these can be used 
in place of the primordial power amplitudes in the parameter set. 

The arena in which CMB theory battles observation is the anisotropy power 
spectrum in multipole space, Fig. 0, which shows how primary Cis vary with 
some of these cosmic parameters. Here Ce = £(t + l)(|(AT/T)£ m | 2 )/(27r). 
Cosmological radiative transfer in the linear regime and therefore the Ci pre- 
dictions as the cosmological parameters are varied is on a firm theoretical 
footing; for reviews, see El [jjl 0]. The Ct sets shown correspond to those in 
Fig. ^. They also are chosen to have the fixed cosmological age of 13 Gyr, 
similar to that inferred from globular cluster ages. Bond and Jaffe |?J used all 
of the current CMB data to test such sets for tilted ACDM, hACDM, oCDM 
sequences, with and without gravity wave contributions and wide variation in 
ft A, Ho, n s , etc. and for ages from 11 to 15 Gyr. Results quoted here on 
parameters from current CMB data and combining CMB and LSS data are 
from that work. 

The C/s are normalized to the 4-year ci?7ir(53-|-90-)-31)(A-|-B) data, which 
fixes it to within about 7% in amplitude. DMR is fundamental to analyses of 
the VLSS region and ULSS region, and is the data set that is the most robust 
at the current time. Even with the much higher precision MAP and Planck 
experiments we do not expect to improve the results on the COBE angular 
scales greatly because the 4-year COBE data has sufficiently large signal-to- 
noise that one is almost in the cosmic variance error limit (due to realization 
to realization fluctuations of skies in the theories) which cannot be improved 
upon no matter how low the experimental noise. 

The "beyond" land in Fig. [j] is actually partly accessible because ultra- 
long waves contribute gentle gradients to our CMB observables. The DMR 
data is well suited to testing whether radical broken scale invariance results 
in a huge excess or deficit of power in the COBE fc-space band, e.g., just 
beyond fc _1 ~ Hq 1 , but this has not been much explored. The remarkable 
non-Gaussian structure predicted by stochastic inflation theory would likely 
be too far beyond our horizon for its influence to be felt in the DMR data. The 
bubble boundary in hyperbolic inflation models may be closer and its influence 
shortly after quantum tunneling occurred could possibly have observable con- 
sequences for the CMB. Theorists have also constrained the scale of topology 
in simple models. For compact spatial manifolds (which may have fl < 1 as 
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Figure 3: The Ce anisotropy bandpower data up to summer 1998 compressed 
to 9 bands using the methods of || are compared with various 13 Gyr model 
sequences (left to right): (1) Hq from 50 to 90, Ha, to 0.87, for an untilted 
ACDM sequence; (2) n s from 0.85 to 1.25 for the H a = 70 ACDM model 
(Qa = .66) - dotted is 0.85 with gravity waves, next without, upper dashed 
is 1.25, showing visually why n s is found to be unity to within 5%; (3) fish 2 
from 0.003 to 0.05 for the H = 50 sCDM model; (4) H Q =50 sequence with 
neutrino fractions varying from 0.1 to 0.95; (5) an isocurvature CDM sequence 
with positive isocurvature tilts ranging from to 0.8; (6) ilsh 2 from 0.003 to 
0.05 for the H Q = 70 ACDM model; (7) H from 50 to 65, fl k from to 0.84 for 
the untilted oCDM sequence, showing the strong ^-shift of the acoustic peaks 
with 17fc ; (9) sample defect Cis for textures, etc. from Jl5| - cosmic string Ce's 
from [ [l6| are similar and also do not fare well compared with the current data. 
The bottom right panel is extended to low values to show the magnitude of 
secondary fluctuations from the thermal SZ effect for the ACDM model. The 
kinematic SZ Ce is significantly lower. Dusty emission from early galaxies may 
lead to high signals, but the power is concentrated at higher £, with a weak tail 
because galaxies are correlated extending into the I < 2000 regime. Forecasts 
of how accurate Ce will be determined for an sCDM model from MAP (error 
bars growing above I ~ 700) and Planck (small errors down the Ce damping 
tail) are also shown. 



well as > 1), the wavenumbers have an initially discrete spectrum, and are 
missing ultralong waves, limited by the size of the manifold. Usually isotropy 
is radically broken in such models, resulting in CMB pattern formation which 
allows much stronger constraints to be imposed than results just from the lack 
of ULSS power; e.g., Q find for flat equal-sided 3-tori, the inscribed radius 
must exceed l.l(2H^ r f= 6600h _1 Mpc from DMR at the 95% confidence 
limit, with the weaker > 0.7(2H ~ 1 ) constraint for asymmetric 1-tori. It is 
also not as strong if the platonic-solid-like manifolds of compact hyperbolic 
topologies are considered, though the overall size of the manifold should be at 
least of order the last scattering surface radius |fj"l[ . 

The primoridial spectral tilt n s is surprisingly well determined using CMB 
data alone, as the upper right panel of Fig. [}] suggests. If as is marginalized for 
the tilted ACDM sequence with Hq=50, with DMR only the primordial index 

I HQ | no 

is n s = 1.02_ 25 with no gravity waves, vt=0, and 1.02_; 1 | with gravity waves 
and v t =i> s , rather encouraging for the nearly scale invariant models preferred 
by inflation theory. For the 13 Gyr tilted ACDM sequence, when all of the 
current CMB data are used get 1.02t;^ for H = 50 (and fi A = 0, the 
tilted sCDM model sequence) and l.OOj/^ for H = 70 (and fl A = 0.66). 
Marginalizing over Hq gives 1.01+'q| with gravity waves included, 0.98j?Qg if 
they are not. The marginalized 13 Gyr tilted oCDM sequence gives 1.00j/Qg. 

Constraints on such "global parameters" as average curvature, Ho and SIa 
from COBE data alone arc not very good, and the situation for the ACDM 
sequences is not that much improved with all of the CMB data. After marginal- 
izing over all n s , get Hq < 75 at Ict, but effectively no constraint at 2a. 
The strong dependence of the position of the acoustic peaks on fife means that 
the oCDM sequence is better restricted: fttot ~ -7 is preferred; for the 13 Gyr 
sequence this gives Hq w 53 and for the 11 Gyr sequence Ho « 65. 

1.4 Cosmic Parameters and LSS (+CMB) 

We have always combined CMB and LSS data in our quest for viable models. 
DMR normalization determines eg to within 7%, and comparing with the <r 8 ~ 
0.6f2~°' 56 target value derived from cluster abundance observations severely 
constrains the cosmological parameters defining the models. In Fig. ||, this 
means the COBE- normalized V$(k) must thread the "eye of the needle" in 
the cluster-band. 

Similar constrictions arise from galaxy-galaxy and cluster-cluster clustering 
observations: the shape of the linear V p must match the shape reconstructed 
from the data. The one plotted in the panels of Fig. || is from ||]. The 
clustering observations are roughly compatible with an allowed range 0.15 < 
r + v s /2 < 0.3, where the oft-used T = T eq e -(^(i+^(2h) 1/2 )-o.06) ; T ^ w 
fl nr h [O er ./(1.68r2 7 )] _1 / 2 , characterizes the density transfer function shape. 
The reason T eq looms so large is that it parameterizes the scale of the horizon 
when the energy density in nonrelativistic and relativistic matter are the same, 



k H l q = 5r- 1 h _1 Mpc. The sCDM model with fl B = 0.03 has T m 0.5. The 
phcnomcnological T/T eq correction factor cannot fully parameterize the effects 
of increasing f2s on the power spectrum, which has the "Sakharov oscillations" 
evident in Fig || in the ABCDM models. 

To get r + zAj/2 in the observed range one can: lower h, lower Q„ r (ACDM, 
oCDM), raise O er , the density parameter in rclativistic particles (1.68f2 7 with 
3 species of massless neutrinos and the photons), e.g., as in rCDM, with a 
decaying v of lifetime r d and V rj 1.08O nr h(l + 0.96(m I ,r d /keV yr) 2 / 3 )" 1 ^. 
raise fi_e; tilt v s < 0, e.g., for sCDM parameters 0.3 < n s < 0.7 would 
be required. Adding a hot dark matter component gives a power spectrum 
characterized by more than just T since the neutrino damping scale enters as 
well as kiieq, as is clear from Fig. |[ In the post-COBE era, all of these models 
that lower T + v s /2 have been under intense investigation to see which survive 
as the data improves, if any. 

The deviation in the shape derived using the APM data 0] from the simple 
r law form may require more baroque models, such as tilted hACDM models 
(see e.g., 0,0). 

Combining LSS and all CMB data gives more powerful discrimination 
among the theories. The approach used to add LSS information to the 
CMB likelihood functions was to design prior probabilities for V + v s j2 and 
i78^nr 6 j reflecting the current observations, but with flexible and generous non- 
Gaussian and asymmetric forms to ensure the priors can encompass possible 
systematic problems in the LSS data. For example, their choice for trgf2^ 6 
was relatively flat over the 0.45 to 0.65 range. 

Using all of the current CMB data and the LSS priors, for the 13 Gyr 
ACDM sequence with gravity waves included, Q get n s = 1.00±;o3 and H = 
72 ± 3 (Oa ~ 0.7), respectively, when Hq and n s are marginalized; with no 
gravity waves, 0.96^ Qg and Hq = 70 ± 3 are obtained; and for an hACDM 
sequence, with a fixed ratio Qhdm/^nr = 0.2 for two degenerate massive neu- 
trino species, n s « 0.97^'q 2 , an d Hq « 57^ are obtained, revealing a slight 
preference for £1a ~ 0.3. A slight increase in age above 13 Gyr lowers Ho to 
the perhaps more preferable 65. 

For the 13 Gyr oCDM sequence, best fit CMB-only models have as too low 
compared with the cluster abundance requirements, so the joint CMB+LSS 
maximum likelihood is substantially below that for ACDM, and severely chal- 
lenged by the data just as it is for the SNIa data. 

The isocurvature CDM models with tilt > 0, e.g., |lj, |l3|], can more 
or less fit the V p shape data in Fig. 0, but are certainly challenged by the 
current CMB data in Fig. |. ABCDM models with an increasing value of Qb 
do have problems with the CMB data, but it can be partly overcome with tilt; 
the deviations in the V p shape from the observed are rather too pronounced 
though {e.g., 10). Calculations of defect models {e.g. strings and textures) 
give CiS that do not have the prominent peak that the data seem to indicate 
[ fL5[ [HI , as Fig. |^ shows. These are only a few of the many examples in 
which the CMB+LSS data has already narrowed our attention enormously in 



structure formation model space. 



2 Nonlinear Probes 

We have mentioned that the magnitude of the linear <J p {k) is a good monitor of 
the sort of dynamics at the resolution that the /c-scale defines. For a p (k) < 0.1 
we are solidly in the linear regime. The weakly nonlinear regime between 0.1 
and cr p ~ 0.7 defines the fc-band largely responsible for the cosmic web, while 
the regime between 0.7 and 2 encompasses most of the virialized objects. 

Over time there have been a number of approximation techniques proposed 
to mimic LSS dynamics. They usually involve the weakly nonlinear regime 
and were designed to be significantly faster than the full A^-body calculations. 
However, with the advent of very fast computing, unless they had some sort 
of semi-analytic counterpart to enhance the depth of understanding such tech- 
niques are of ever-diminishing importance. The most venerable and still useful 
is the Zeldovich approximation, i.e., first order Lagrangian perturbation the- 
ory. Second order Lagrangian perturbation theory is much more accurate, 
but also harder to implement. The Zeldovich approximation with sticking at 
caustic formation is called adhesion theory p8[ , utilizing Burgers equation for 
the velocity. It successfully describes the architecture of voids, filaments and 
clusters but does not deal with interior dynamics of "stuck structure" and 
has not been very useful for mass estimates. The frozen potential and frozen 
flow approximations are attempts to avoid full potential calculations, but are 
instrinsically numerical, and do not work that well. Sec |lS| for a review. 

What has been applied with great success is weakly nonlinear perturba- 
tion theory, both Eulerian and Lagrangian, much of it by our French hosts 
(see e.g., p0|). Some of the milestones in this effort are: accurate prediction 
of the one-point distribution function (PDF) for the overdensity S at "tree 
level" c.f. N-body results (2l], |22j; analytic expressions for the cumulants 
S n = (S n )/(°±N L ) n ~ 1 t0 " tree ordei "" s n ~ 0(a 2n ~ 2 ) agree well with TV- 
body results p3[ ; calculation of the variance (5 2 ) = a 2 NL to "one-loop" order 
[ p4| p5{ and, in a tour de force, S3 j26j and the bispectrum |^7| to this order 
(the az NL order correction to the tree level results). 

As one heads into somewhat stronger nonlinearity, a p > 0.7, analytic meth- 
ods to deal with object collapse are heavily utilized. The hugely popular, 
trivial-to- implement, Press-Schechter (1976) method, as modified by |35| to 
the excursion set theory: at each point in space, 5h{k) is allowed to random 
walk as the resolution a p increases, i.e., as the scale decreases. Once the 
5i(fc) reaches a redshift-dependent threshold f c (z), the "absorbing barrier", 
that point is said to have collapsed at that redshift and is assigned to a halo 
with mass M corresponding to the scale k. The flaw is evident: nearby 
points that belong to the same halo will be assigned to different mass objects. 
Thus the amazement, and delight, in the community that the one-point distri- 
bution, n(M)d\nM, so derived fits that of A^-body group catalogues so well 



[||[ |36[ even below M*. So do such interesting constrained mass func- 
tions as n(M2, Z2\M\, Zx) giving the number density of objects of mass M% at 
redshift Z2 given that the region is already within one of mass M\ at redshift 

Zx. 

The peak-patch picture p4j described below is the natural generalization 
of the Press-Schechter method to include non-local effects, spatial correlations 
and more natural ways of assigning mass to halos. It is also the natural 
generalization of BBKS single-filter peaks theory j3?j to allow a mass spectrum 
and solve the cloud-in-cloud (i.e., peak-in-peak) problem. In it, the threshold 
f c becomes a function of the tidal environment the peak patch finds itself in. 
It is also at the heart of our cosmic web picture. 

Another approach to the nonlinear evolution of V p (k) is that of Hamilton et 
al. [p8| , where a Lagrangian wavenumber &l of linear theory is mapped onto a 
nonlinear wavenumber k^L through k^ L w (I+^atl) -1 ^^ 1 , where 5nl is the 
average nonlinear overdensity in a mean field region of scale fc^L' J us ^ wna t y° u 
would expect if you took a spherical Lagrangian radius and compressed it to an 
Eulerian radius while conserving the mass. A mapping from the linear power 
spectrum V p l to the nonlinear is defined by Ppjvi(fcjvi) = fn^f pi (fee)), where 
f n denotes a suitably fit nonlinear function, which, however, is dependent upon 
the shape of the linear power spectrum |^9) . The VIRGO consortium tested 
the fn fit given by p0[ with their A-body runs, and find accurate mapping 
from V p l to power levels V p nl up to ~ 1000. It was used by B to take 
the nonlinear galaxy power spectrum derived e.g., from the APM survey and 
make an estimate of the linear one: this is what we adopted for comparing 
with theory in Fig. ^. 

There is of course the direct numerical approach to nonlinear physics. We 
describe the various A-body and hydro methods currently employed in cos- 
mology. We then turn our attention to the peak-patch picture and the cosmic 
web. 

2.1 N-Body and Gas Simulations 

The ITP Cluster Comparison of Cosmological Hydro Codes was a "homework 
problem" assigned to simulators at the extended workshop on Galaxy Forma- 
tion and Cosmic Radiation Backgrounds held at the Institute for Theoretical 
Physics in Santa Barbara in 1995. Calculations were finished in 1996/97 and 
the paper was submitted in 1998 |3]}| . A constrained single peak field was 
used for initial conditions to ensure that a massive COMA-like cluster would 
form at the centre of the simulation. The results provide a good snapshot 
of the A-body and hydro methods currently in use in cosmology, although 
improvements in all of the codes have been made since the test. The people 
participating and the codes used are listed in Table 

In spite of the huge variation in code types, CPU hours devoted, memory 
required, mass resolution and spatial resolution used, there was surprisingly 
good agreement among the calculations: excellent in the densities of gas and 



ITP Cluster Comparison of Cosmological Hydro+N-body Codes |3l[] 



Group 


Method 


CPU 


machine 


storage 


Bond & Wadsley 


SPH+P 3 MG 


119hr 


DECa 


100MB 


Bryan & Norman 


PPM+PM 


200 


SGI PowCh 


500 


Cen 


TVD+PM 


5312 


IBM Sp2 


4400 


Evrard 


SPH+P 3 M 


320 


HP375 


17 


Gnedin 


SLH+P 3 M 


136 


SGI PowCh 


90 


Jenkins, Thomas & Pearce 


SPH+AP 3 M 


5000 


Cray-T3D 


512 


Owen & Villumsen 


ASPH+PM 


40 


Cray-YMP 


106 


Navarro 


SPH+Direct 


120 


SparclO+Grape 


75 


Pen 


MMH+MMPM 


480 


SGI PowCh 


900 


Steinmetz 


SPH+Direct 


28 


SparclO+Grape 


22 


Couchman 


SPH+AP 3 M 


77 


DECa 


95 


Yepes & Klypin 


FCT+PM 


350 


Cray-YMP 


480 


Warren & Zurek 


Tree 


15360 


Intel-A 


1000 



GRAVITY SOLVERS: Direct refers to the brute force direct sum over all pairs to 
get the force. Tree refers to organizing distant particles into groups using a tree. 
PM refers to particle-mesh, using an FFT to compute the gravitational potential. 
MMPM is PM with a Lagrangian mesh distorted to move with the flow. P 3 M refers 
to particle-particle + PM, using FFTs to get the large scale force, and direct sum to 
get the small range force. AP 3 M is P 3 M with adaptive mesh refinement, so not as 
much costly P-P is needed. The preferred method of Bond & Wadsley is treeP 3 M, 
PM |^| for long-range forces with a tree-like P-P for the short-range. It is twice 
faster than the SPH-P 3 MG code used in our cluster comparison, with more accurate 
forces and half the storage requirement. (P 3 MG refers to particle-particle plus a 
multigrid solution to the Poisson equation to get the long range force.) Warren and 
Zurek did DM only, but a very large massively parallel benchmark calculation. 
LAGRANGIAN HYDRO SOLVERS: These methods have resolution following the 
flow, hence increasing in high density regions. SPH is Monaghan's smooth particle 
hydrodynamics with spatially variable resolution; ASPH is aspherical SPH, using 
anisotropic kernels to treat asymmetric collapses better; MMH is moving-mesh hy- 
dro using a pseudo-Lagrangian grid-based technique; SLH is grid-based softened La- 
grangian hydro. The Jenkins et al. simulation was a highly parallelized calculation 
using the SPH code Hydra; Couchman used it on a serial machine. 
EULERIAN GRID-BASED HYDRO SOLVERS: PPM is the piecewise parabolic 
method; TVD refers to the Total Variation Diminishing hydro scheme; FCT refers to 
the Flux-Corrected- Transport scheme. For these Eulerian schemes to ultimately be 
effective in such compressible media as we deal with in cosmology, AMR (adaptive 
mesh refinement) must be used, which the PPM code has (an early single AMR 
version was used). 

dark matter; good in the gas temperature, entropy, pressure, X-ray luminosity 
density, and gas to dark matter ratio; not as good in the highly dynamical 
earlier phases before cluster virialization. 

The largest N-body simulation to date has been done by the VIRGO consor- 
tium, 1 billion particles using a P 3 M code of a periodic region (2000 h _1 Mpc) 3 
in size, on a CRAY T3E with 512 CPUs, of a V = 0.21, cr 8 = 0.6 rCDM model. 



The mass resolution denned by the mass per particle was 4 x 10 12 Mq and the 
gravitational softening was 200 kpc, chosen to ensure a do-able calculation by 
avoiding a large P-P slowdown. The texture of the cosmic web is very evident 
in the results ffl . 

We have adopted a different approach we refer to as importance sampling 
[ [32"| , in which we do a number of simulations of constrained field initial con- 
ditions, chosen to sample most efficiently the statistical distribution of what- 
ever quantities we wish to probe. It is as if we do a big-box simulation, and 
go around sampling this patch and that until we get statistical convergence. 
From this point of view, the big box simulations oversamplc the patches with 
large scale fluctuations near the rms, yet may undersample the rare protosu- 
percluster or protosupervoid regions. If the periodic simulation volumes are 
not big enough so that there is negligible density power on the scale of the 
fundamental mode, k^ 1 = L/(2ir), then the big-box results will be in error. 
This is not a problem for the VIRGO simulation, but is a major headache for 
simulations of the Lyman a forest where cr p (k) changes slowly, and even more 
so for simulations of the "1st *" region. We design our simulation patches to 
avoid this For many purposes, we have found it adequate to calculate 

"shearing patches" , constrained by specifying the tidal tensor at the origin 
smoothed on a few scales. However, we also may be interested in simulating 
complex regions, where more controlling constraints are used, e.g., associated 
with many cluster-scale density peaks. Our example of a dense supercluster is 
of this form. 

We now describe some aspects of simulation design. We first decide on 
the mass resolution we wish to achieve. This is set by the lattice spacing 
of particles on the initial high resolution grid, ax , chosen to be ~ 2 Mpc 
to ensure that the waves needed to treat the target objects forming in the 
simulation will be adequate, in this case clusters. Next we need to determine 
the spatial resolution, of the gas and of the gravitational forces, preferably 
highly linked. In Lagrangian codes like we use, this varies considerably, being 
very high in cluster cores, moderate in filaments, and not that good in voids. 
Here we wanted to get good calculations of the cluster cores, so we wanted 
our resolution to be in excess of ~ 40 kpc, and the value we get is about 
30 kpc. Given the target resolution, we then have to determine how large 
the high resolution part of the simulation volume is by CPU limitations on 
the number of particles we can run in the desired time. This may mean the 
high resolution volume may distort considerably during the simulation. To 
combat this, layers are added of progressively lower resolution volumes, to 
ensure accurate large scale tides/shearing fields operate on the high resolution 
patch. For the calculations shown, the High Resolution region of interest (grid 
spacing az,, 50 3 sphere) sits within a Medium Resolution region (2az,, 40 3 ), in 
turn within a Low Resolution region (2az,, 32 3 ). The mean external tide of 
the entire patch is linearly evolved and applied during the calculation. 

Thus our sample simulation has a high resolution 104 Mpc diameter patch 
with 50 3 /2 gas and 50 3 /2 dark matter particles (initial grid spacing of 100 



comoving kpc), surrounded by gas and dark particles with 8 times the mass 
to 166 Mpc, in turn surrounded by "tidal" particles with 64 times the mass to 
266 Mpc. Bigger calculations are easily do-able without resorting to massive 
parallelization; e.g., a 100 3 /2 gas and 100 3 /2 dark matter particle calculation 
with a total of 1.6 million particles including the medium and low resolution 
regions takes about two weeks of SGI Origin 2000 single processor time. 

The mass resolution limits the high k power of the waves that can be laid 
down in the initial conditions (Nyquist frequency, it/cll), but for aperiodic 
patches there is no constraint at the low k end: we use the FFT for high 
k, but a power law sampling for medium k, and a log k sampling for low k, 
the latter two done using SlowFT, i.e., a direct sum over optimally-sampled 
k values, with the shift from one type of sampling to another determined 
by which gives the minimum volume per mode in fc-space. By contrast, the 
periodic big box calculations are limited by the fundamental mode. 

The numerical method we adopt is the cosmological SPH+treeP 3 M code of 
Wadsley and Bond Jij]| . The treeP 3 M technique is a fast, flexible method for 
solving gravity and can accurately treat free boundary conditions. The code 
includes photoionization as well as shock heating and cooling with abundances 
in chemical (but not thermal) equilibrium, incorporating all radiative and 
collisional processes. The species we consider are H, H + , He, He + , He ++ and 
e _ . The simulations were run from z = 30 to z = 0. Parameters for the 
sCDM, ACDM and oCDM cosmologies we have run are given in Fig. |]. 

Our a,L resolution is 2 times better than the VIRGO simulation, and our 
best spatial resolution is 6 times better. On the other hand, we are simulating 
only 10~ 5 of the VIRGO volume at high res: even with clever sampling, being 
computer-cheap does have its drawbacks; though this is a rare event, about a 
thousand similarly compact superclusters would be within the VIRGO volume. 
We are now running simulations at twice the resolution to ensure accuracy for 
X-ray emission in the filaments, not just in the clusters. Our resolution is 
adequate for making predictions for signals from the supercluster complex in 
the Sunyaev-Zeldovich effect (Fig. ||) and gravitational lensing (described in 
our companion paper [p9| ). 

2.2 The Peak-Patch Picture and the Cosmic Web 

We approach the connected ideas we loosely call Cosmic Web theory (Bond, 
Kofman and Pogosyan 1996, BKP) and the Peak-Patch Picture (Bond and 
Myers 1996, BM, where detailed references to the pre-1995 work quoted here 
are given) through a historical path that includes an outline of the relevant 
terminology. In 1965, Lin, Mestel & Shu showed that a cold triaxial collapse 
implied an oblate "pancake" would form. In 1970, Zeldovich developed his 
famous approximation and argued that pancakes would be the first structures 
that would form in the adiabatic baryon-dominated universes popular at that 
time. Generally for a cold medium, there is a full non-linear map: x(r, t) = r — 
s(r,i), from Lagrangian (initial state) space, r, to Eulerian (final state) space, 
x. The map becomes multivalued as nonlinearity develops in the medium. 



Figure 4: COMPACT SUPERCLUSTERS: Dark matter density in a patch 
~ 100 Mpc (comoving diameter) across at redshifts 2=1, 0.5 and (from top 
to bottom) of three constrained-field supercluster simulations with differing 
cosmologies, with HR region 104 Mpc, MR region 166 Mpc and LR region 
266 Mpc, with wave coverage to k^ 1 = 1000 Mpc, the very long waves enter- 
ing through a self consistent mean tidal field. The cosmologies shown from 
left to right are standard CDM (er 8 =0.67, h=0.5), ACDM (cr 8 =0.91, h=0.7, 
a ir =0.3349) and an Open CDM model (ct 8 =0.91, h=0.7, fi„ r =0.3689). The 
initial conditions were constructed using constraints from the peak patches ei- 
ther within the region or those exerting a significant tidal force on it. e.g., 17 
peak patches (found at z = in a 400 Mpc simulation) were used for sCDM. 
The constraints imply the simulations have similar, though not identical, pat- 
terns. Whereas the middle panel ACDM simulation had 29, 10 and 2 clusters 
with mass above 3 x 10 14 M Q at redshift 0, 0.5 and 1 and oCDM had 20, 12 
and 5, sCDM had 13, 3 and none, a dramatic statement of the very different 
redshift evolution of the cluster system for sCDM. 
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Figure 5: SUNYAEV-ZELDOVICH MAP for the supercluster region seen at 
redshifts 1, 0.5 and 0.2 subtending the angles shown. The observing wavelength 
was taken to be in the Rayleigh- Jeans region of the spectrum (so AT/T = — 2y 
here). The core red regions are above 32 x 10 -6 , and the dark contours 
surrounding the white are at 2 x 10~ 6 . These levels are now accessible to 
ground based instrumentation. Blank field SZ surveys using interferometers 
and bolometer experiments promise to revolutionize our approach to the clus- 
ter system, especially at z > 0.5. Although the contours do go into the far 
field of the cluster region because they probe the electron pressure, one needs 
to go below 10 -6 to fully detect the filaments. The effect is more pronounced 
in the ACDM and oCDM universes than in the sCDM model. The analogue 
of this figure in weak lensing is shown in p9[ . 



It is conceptually useful to split the displacement field, s = Sb + s/, into a 
smooth quasilinear long wavelength piece s& and a residual highly nonlinear 
fluctuating field Sf. As we have discussed, if the rms density fluctuations 
smoothed on scale Rb, a p (Rb), are < 0(1/2), the s&-map is one-to-one (single- 
stream) except at the rarest high density spots. In the peak-patch approach, 
Rb is adaptive, allowing for dynamically hot regions like protoclusters to have 
large smoothing and cool regions like voids to have small smoothing. If D(t) 
is the linear growth factor, then Sb = D(t)sb(r,0) describes Lagrangian linear 
perturbation theory, i.e. the Zeldovich approximation. The large scale peculiar 
velocity is ~Vpb — —a(t)sb(r, t). What is important for us is the strain field (or 
deformation tensor): 

1 ds ds 

e b .rj{r) = -3(7^ + = -J2 XvAfl ' vAfl " A ' where 

At,3 = -^r(! + 3e « +Pv), K2 = -^(1 -2p v ), Xvi = ^p(l + 3e„ -p v ), 

and S^b = —e\ i is the smoothed linear overdensity, which we often express 
in terms of the height relative to the rms fluctuation level a p (Rb), Vb = 
8 Lb I 'a p(Rb) ■ The deformation eigenvalues are ordered according to A„3 > 
^v2 > A-ui and n v A denote the unit vectors of the principal axes. In that 
system, xa — ta(1 — \ v A(r,t)) describes the local evolution.^ 

The strain tensor is related to the peculiar linear tidal tensor by g§r^ = 
—4irGp nr a 2 eb t ij, where $ is the peculiar gravitational potential, and to the 
linear shear tensor by the time derivative e^-. The anisotropic part of the 
shear tensor has two independent parameters, the ellipticity e v (always posi- 
tive) and the prolaticity p v . 

Doroshkevich (1973) and later Doroshkevich & Shandarin (1978) were 
among the first to apply the statistics of Gaussian random fields to cosmol- 
ogy, in particular of X v a, at random points in the medium (see also p9j). 
Arnold, Shandarin and Zeldovich (1982) made the important step of applying 
the catastrophe theory of caustics to structure formation. This work sug- 
gested the following formation sequence: pancakes first, followed by filaments 
and then clusters. This should be compared to the BKP Web picture forma- 
tion sequence: clusters first, followed by filaments and then walls. BKP also 
showed that filaments are really ribbons, walls are webbing between filaments 
in cluster complexes, and that walls are not really classical pancakes. For 
the Universe at z ~ 3, massive galaxies play the role of clusters, and for the 
Universe at z ~ 5 more modest dwarf galaxies take on that role. 

The Web story relies heavily upon the theory of Gaussian random fields as 
applied to the rare "events" in the medium, e.g., high density peaks. Salient 
steps in this development begin with BBKS (1986) J37J, where the statistics 

2 The overdensity (1 + 6 z )(r, t)= |(1 - D(t)A„ 3 )(l - D(t)X v2 )(l - DftJA^)!" 1 in a Zel- 
dovich map explodes when the largest eigenvalue D(t)X v s reaches unity (fold caustic forma- 
tion); a pancake develops along the surface n„3 • V r A„3 = 0. In hierarchical models, classic 
Zeldovich pancakes are not relevant for structure formation. 



of peaks were applied to clusters and galaxies, e.g., the calculation of the 
peak-peak correlation function, £, p k, P k- In a series of papers, Bond (1986-90) 
and Bond & Myers (1990-93) developed the theory so that it could calculate 
the mass function, n(M)d\nM. It was also applied to the study of how shear 
affects cluster alignments (e.g., the Binggeli effect), and to Lya clouds, 'Great 
Attractors', giant 'cluster-patches', galaxy, group and cluster distributions, 
dusty PGs, CMB maps and quasars. This culminated in the BM "Peak-Patch 
Picture of Cosmic Catalogues" . 

We briefly describe the BM peak-patch method and how it is applied to 
initial conditions for simulations such as in Fig. [|. We identify candidate 
peak points using a hierarchy of smoothing operations on the linear density 
field Sl- To determine patch size and mass we use an ellipsoid model for the 
internal patch dynamics, which are very sensitive to the external tidal field. 
A byproduct is the internal (binding) energy of the patch and the orientation 
of the principal axes of the tidal tensor. We apply an exclusion algorithm 
to prevent peak-patch overlap. For the external dynamics of the patch, we 
use a Zeldovich-map with a locally adaptive filter (R p k) to find the velocity 
Vpfc (with quadratic perturbation theory corrections sometimes needed) . The 
peaks are rank-ordered by mass (or internal energy). Thus, for any given 
region, we have a list of the most important peaks. By using the negative of 
the density field, we can also get void-patches. 

The peak-patch method allows efficient Monte Carlo constructions of 3D 
catalogues; gives very good agreement with ./V-body groups; has an accurate 
analytic theory with which to estimate peak properties, (e.g., mass and bind- 
ing energy from mean-profiles, using f c (e v ), (e v \v p k))\ and handles merging, 
with high redshift peaks being absorbed into low redshift ones. 

BKP concentrated on the impact the peak-patches would have on their 
environment and how this can be used to understand the web. They showed 
that the final-state filament-dominated web is present in the initial conditions 
in the &Lb pattern, a pattern largely determined by the position and primordial 
tidal fields of rare events. BKP also showed how 2-point rare-event constraints 
define filament sizes (see Fig. |J). The strongest filaments are between close 
peaks whose tidal tensors are nearly aligned, a binary molecule image with 
oriented peak-patches as the atoms. Strong filaments extend only over a few 
Lagrangian radii of the peaks they connect. They are so visually impressive 
in Eulerian space because the peaks have collapsed by about 5 in radius, 
leaving the long bridge between them, whose transverse dimensions have also 
decreased. The reason for the strong filaments between aligned peaks is that 
the high degree of constructive interference of the density waves required to 
make the rare peak-patches, and to preferentially orient them along the 1-axis, 
leads to a slower decoherence along the 1-axis than along the others, and thus 
a higher density. Membrane walls are the less dense webbing between the 
filaments, a 3,4, . . .-peak molecular image, also shown in Fig. ^. And void- 
patches are the inverse of peak-patches in the initial conditions, but mapped 
by nonlinear dynamics into the dominant volume elements. 




Figure 6: These plots illustrate the molecular picture of large scale structure, 
with "bonds" bridging clusters. The initial conditions have been smoothed and 
Zeldovich-mapped. (Full iV-body maps look similar.) Upper panels show a 
two-point mean field (5i|2pks) constrained by two oriented clusters separated 
by 40h _1 Mpc, fully aligned and partially aligned. The next 4 panels show 
three-point and four-point mean fields for different peak-patch orientations 
taken from a simulation. Notice the lower density contrast webbing between 
the filaments. One of the supercluster simulations had initial conditions con- 
strained by the clusters in this region. 



The Cosmic Web and Peak-Patch pictures provide powerful language for 
understanding the structure and evolution of not only cluster systems, but 
also galaxy formation and Lyman alpha absorption systems at high redshift. 
The web theory predicts the basic structural components of the IGM as a 
function of scale, epoch and cosmology. For density contrasts 6 >100, the 
rare-events at z ~ 3 are massive galaxies, observed as damped Lya absorbers, 
while the typical collapsed objects are dwarf galaxies responsible for Lyman 
Limit and metal line systems. The medium is visually dominated by 5 ~ 
5 — 10 filaments, bridging massive galaxy peaks, with smaller scale filaments 
within the larger scale ones bridging smaller dwarf galaxies contributing the 
most to the Nhi ^3 10 Lya forest. Peak-patch galaxy catalogues with 
all large scale waves included covering very large volumes of space can be 
constructed relatively cheaply for comparison with observations GJ, e.g., with 
the impressively grand Steidel et al. structures at z ~ 3.1 Q, a task beyond 
current A-body computational capability. A natural byproduct, of course, is 
the increasing bias of halos with increasing velocity dispersion or mass. 

We have used field realizations constrained to have interesting multiple 
peak/ void structures as initial conditions for high resolution numerical sim- 
ulations, e.g. the strong filament of galaxies in |32|. In this paper, we had 
in mind mimicking the remarkable Shapley concentration of rich clusters with 
the simulations shown in Fig. ^. All peak patches tidally influencing a region 
that had one of the highest smoothed densities of protocluster-patches in ini- 
tial conditions for a (400 Mpc) 3 A-body simulation were used to set up new 
initial conditions with the same basic webbing pattern for our hydro calcula- 
tions. For the sCDM case, for example, 17 peaks were used, including the four 
nearby clusters shown in the lower panel of Fig. ^| which exhibited such strong 
filamentary bridging. These reside in the core region of Fig. ^. 

We have descended from ethereal realms of early universe particle physics, 
through the calm of linear CMB and matter transport, into the elegance of the 
weakly and even strongly nonlinear collisionless gravitational problems; and 
finally into gastrophysics, where the "subgrid physics" of star formation, mul- 
tiphase ISM, feedback from winds, explosions, ionization fronts, etc. must be 
injected into the simulatable scales, providing work for decades to come. Even 
so, with the amazing databases soon to appear and our developing theoretical 
tools, we might just possibly answer how large scale structure formed. 
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